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We study real-space condensation phenomena in a type of classical stochastic processes (site- 
particle system), such as zero-range processes and urn models. We here study a stochastic process 
in the Ehrenfest class, i.e., particles in a site are distinguishable. In terms of the statistical mechan- 
ical analogue, the Ehrenfest class obeys the Maxwell-Boltzmann statistics. We analytically clarify 
conditions for condensation phenomena in disordered cases in the Ehrenfest class. In addition, we 
discuss the preferential urn model as an example of the disordered urn model. It becomes clear 
that the quenched disorder property plays an important role in the occurrence of the condensation 
phenomenon in the preferential urn model. It is revealed that the preferential urn model shows 
three types of condensation depending on the disorder parameters. 

PACS numbers: 75.10.Nr, 02.50.Ey, 64.60.-i, 05.40.-a 



I. INTRODUCTION 

Condensation phenomena have been widely observed 
in various stochastic processes and physical contexts, 
such as mass transport models P, [1, Q j traffic flow Q , 
condensation of edges in complex networks [s] . It is no- 
table that some of the condensation phenomena occur 
in nonequilibrium systems, such as asymmetric simple 
exclusion process (ASEP) and zero-range process (ZRP) 
. In general, it is difficult to study such nonequilibrium 
systems, but the steady state properties in a class of the 
stochastic processes can be treated analytically because 
their partition functions are written in simple factorized 
forms. Such a simple form of the partition function en- 
ables us to study the condensation phenomenon analyti- 
cally. 

Mass transport models, such as the ZRP and urn mod- 
els 0, i, i, [IS [O, [H, [11 , have been widely used in or- 
der to research characteristic features of the condensation 
phenomena. When the condensation arises, a single site 
contains a macroscopic number of particles. Both the 
ZRP and urn models consist of many sites containing an 
integer number of particles. Inter-particle interactions 
are modeled by allowing the particle transition rate to 
depend on the number of particles in each site. Further- 
more, in order to discuss general situations, it is possible 
to consider the case in which the transition rate depends 
also on a site-dependent disorder parameter. The math- 
ematical structures for the ZRP and the urn models are 
almost the same, so that one can treat these models in a 
similar manner. 

In the urn models, there are two types of statistics 
[T^ l : one is the Ehrenfest class in which particles in a site 
are distinguishable; the other is the Monkey class which 
consists of indistinguishable particles. It has been re- 
vealed that the Ehrenfest class corresponds to the statis- 
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tical mechanics with the Maxwell-Boltzmann statistics. 
Monkey class is related to the Bose-Einstein statistics, 
and in general raises real-space condensation phenomena 
1^1 [l^ ■ The real-space condensation phenomenon means 
that a site has a macroscopic number of particles. 

We here note that some of recent researches on com- 
plex networks have dealt with distinguishable compo- 
nents [ii[i3,[il,[il, and a type of condensation phenom- 
ena has been observed numerically in a network model 
in the canonical ensemble, in which the total number of 
nodes and edges do not change [l8| . When we map a com- 
plex network to an urn model, the distinguishable prop- 
erty of edges in the network model corresponds to that of 
particles in the urn model [l^ . We therefore expect that 
the condensation phenomenon is able to occur even in 
the Ehrenfest class (Maxwell-Boltzmann statistics) due 
to the occurrence of the condensation phenomenon in the 
complex networks. Actually, there are a few example of 
the Ehrenfest class, in which the real-space condensation 
phenomenon can arise [l^. [isj. However, in the complex 
networks, quenched disorder in the model is expected to 
be important, and discussion between the quenched dis- 
order and the condensation phenomenon in the Ehrenfest 
class has not been given yet. 

In the present paper, we consider stochastic processes 
with the Maxwell-Boltzmann statistics, i.e., the Ehren- 
fest class, and with the quenched disorder. Conditions 
for the occurrence of the condensation phenomenon are 
clarified. Furthermore, we study one example for the 
condensation phenomenon in the Ehrenfest class with 
the quenched disorder, i.e., the preferential urn model 
proposed in Ref. [loj . We analytically and numerically 
confirm the condensation phenomenon in the preferen- 
tial urn model. It is revealed that there are three types 
of condensation phenomena depending on the disorder 
parameters. 

The construction of the present paper is as follows. In 
Sec. II, we will present examples of the stochastic models 
in the Ehrenfest class (the Maxwell-Boltzmann statistics) 
and explain their partition functions. The conditions for 
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the condensation phenomenon in disordered cases are ex- 
plained in Sec. III. In Sec. IV, we give a concrete example 
for the condensation, i.e., the preferential urn model, and 
discuss three types of the condensation phenomena. We 
also give some numerical results in order to confirm the 
condensation phenomena in Sec. IV. Section V gives con- 
cluding remarks. 



II. MODELS IN THE EHRENFEST CLASS 
A. Examples of models 

We study stochastic processes with a factorized par- 
tition functions. Examples of these stochastic processes 
are the ZRP and the urn models, as described in Sec. I. 
There are two types of statistics 12]: One is the Ehren- 
fest class; the other is the Monkey class. The charac- 
teristics of these two classes are summarized in Table. ID 
The Ehrcnfest class has distinguishable particles, so that 
it corresponds to the Maxwell-Boltzmann statistics. In 
contrast, indistinguishable particles are considered in the 
Monkey class. Here, we give two examples of the stochas- 
tic processes in the Ehrenfest class and with quenched 
disorder. 

• Urn model on a complete graph in Ehrenfest class: 
The urn model on a complete graph consists of L 
sites (urns) and N particles. The density p is de- 
fined by p = N/ L. An disorder parameter hi is 
assigned to each site by a probability density 

The number of particles in site i is denoted by n^. 
When we adopt a Monte Carlo method (the heat- 
bath rule) for the dynamics, a particle is selected 
randomly in one Monte Carlo step, and then trans- 
fered to an urn selected by a probability u' {hi, Ui) 
(Fig. [TJa)) Note that the probability depends 
on the number of particles and the disorder param- 
eter of the arrival site; it does not depend on the 
properties at the departure site. 

• ZRP with the Maxwell-Boltzmann statistics in one- 
dimensional lattice with periodic boundary condi- 
tions: The lattice consists of L sites, and there are 
N particles. Each particle hops out to the nearest 
neighboring site, only in one direction. A particle 
is selected by a probability u{hi,ni) in one Monte 
Carlo step, and hops out (Fig. [Tfb)). 

We note that the most general ZRP has not the 
Maxwell-Boltzmann characteristics. The dynamics of the 
most general ZRP is as follows; "a site is selected by a 
probability u{hi,ni), and a particle in the selected site 
hops out (Fig. [Ijc))." Hence, the particles in a site are 
indistinguishable, so that the most general ZRP corre- 
sponds to the Monkey class. 

In Sec. IV, we will specify an urn model in the Ehren- 
fest class in more detail in order to investigate the 
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FIG. 1: (a) The urn model in the Ehrenfest class (the 
Majcwell-Boltzmann statistics), (b) ZRP in the Ehrenfest 
class, in which each particle can hop with rate u{h,n). (c) 
The most general ZRP (the Monkey class) , in which only one 
particle in a site hops out with rate u{h, n). 



condensation phenomenon in disordered cases with the 
Maxwell-Boltzmann statistics. 



B. Partition function 

For the Monkey class (e.g., the most general ZRP), the 
partition function is written as [l^ 



^Monkey 



m— riL—O 



(1) 



where ^ ^X^i^^i '^i' Kronecker delta, which indi- 
cates that the total number of particles is conserved. 
The weight p{hi, Ui) is calculated from the transition rate 
of the ZRP, u{hi,ni), by p{hi,ni) = 1/ Oj^i ^'(^ii j)- 
For the urn models, we calculate the weight p{hi,ni) as 
p{h,,n,) = lY-U u'{hi,j). 

The partition function of the Ehrenfest class (Maxwell- 
Boltzmann statistics) is a similar one with that of the 
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TABLE I: Comparison between the Ehrenfest class and the Monkey class. 





Ehrenfest class 


Monkey class 


Particles in a site 


distinguishable 


indistinguishable 


Weight of site i 


p{hi,ni)/ni\ 


p{hi,ni) 


Statistical-mechanical analogue 


Maxwell-Boltzmann 


Bose-Einstein 


Condensation phenomenon 


difficult to occur 


easy to occur 



Monkey class; 



E--E 

m— iiL—O 



p{hi,ni) p{hL,nL) 



(2) 



Note that the partition function involves factorials 
{l/(ni!)}, which stems from the distinguishable property 
of particles [l^l- The factor l/(ni!) plays an important 
role for the condensation phenomenon, which will be ex- 
plained in the subsequent section. 



III. CONDENSATION IN DISORDERED CASES 

From the partition function of Eq. Q , the steady-state 
occupation distribution is derived. Since we consider site- 
dependent disordered cases, we analyze it with the replica 
method, as described in Ref. [20]. We here briefly sum- 
marize the analysis pTj . 

Using the partition function of Eq. the probability 
with which site (urn) 1 has k particles is given by 



pation distribution as 



dhi 4'{hi) 
dhi 4>{hi) 



^ I {hx,...JlI,} 

kl \ Zl,n 
P{hi,k) 
k\ 



{/12,...Jil} 



X exp 



(5) 



where (• • • ){/ii./i2,- -,'ii} configurational average, 

and (f){h) a probability distribution for the disordered pa- 
rameters {hi}. The third line in Eq. (O is obtained by 
the use of the self-averaging property 20] . 

In order to calculate (InZ^jv) and (InZ^ ^), we use 
the replica technique. The following replica identity is 
used in order to calculate the configurational average: 



II, ^ 1 = lim 

{h2,...,hL} „j^o 



(6) 



Using the integral representation of Kronecker delta, 



27rj(5(a, b) = f dz 



and defining 



(7) 



Zl 



^ oo oo 



p{hi,ni) p{hL,nL) 



p{hi, k) Z'j^ 



ni! 



N 



kl Zl. 



N 



where we define jy 



= E--E 

n2=0 nL=0 



n2l 



S (^rii,N -kj 



p(/l2,f^2) p{hL,nL) 



(3) 



(4) 



for notational simplicity. We therefore obtain the occu- 



the replicated partition function in the configurational 
average is given by 



dzi f dZr, 



27ri 



X exp |^(L — 1) In ^cxp 



^lni7(/i,z„) 



(8) 



Assuming the replica symmetry (z^ = z) and using the 
saddle-point method, we get the following saddle-point 
equation 



(9) 



and 



(lnZi,Ar){^^_ ~ - (pL + 1) Inzs + \\\B{hx,z^) 

+ {L-l){lnHih,Zs))h + C, (10) 
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where C is a constant. We perform the similar calculation 
for {luZ'j^ j^j), and obtain the same saddle-point equation 
as Eq. ® and 

~ - + 1 - fc)lnzs 

+ {L-l){lnH{h,Zs))h + C, (11) 

where the constants C in (^L,Af){h2,--,'iL} ^^"^ 

{Z']^ Ar){h2 fiL} ^I'G the same. Finally, we get the 

following occupation distribution 

PW = /.M(/.)^^- (12) 

For the condensation phenomenon, the radius of con- 
vergence z in Eq. ([7]) plays an important role. We here 
define the radius of convergence as z — /3. In Eq. ([9]), 
increasing the value of Zg corresponds to the increase of 
the density p{zs). The radius of convergence Zs = P gives 
the maximal value of the density p^'- 

Pc^pl^ /?)) /h{K P)^^ . (13) 

When the density p goes beyond the critical density pc, 
the excess particles condense only on a site. 

It has been revealed that the following conditions are 
needed for the condensation phenomenon [6|: 

• If the radius of convergence f3 is infinite, the con- 
densation phenomena do not occur. 

• If the radius of convergence (3 is finite and pc is 
infinite, one can find a solution of Eq. ([9|) for any 
finite density, so that there is no condensation phe- 
nomena. 

• If both P and pc are finite, one can no longer satisfy 
Eq. © for p > Pc- In this case, we have condensa- 
tion phenomena. 

It is obvious that the condensation phenomenon is dif- 
ficult to occur in the Ehrenfest class, because the fac- 
tor 1/n! makes the summation of Eq. ([7]) converge for 
much choice of p{h,n). Hence, the choice of the transi- 
tion probability (or p(h,n)) would be important for the 
condensation phenomenon. 

IV. SPECIFIC EXAMPLE: DISORDERED 
PREFERENTIAL URN MODEL 

A. Definition of tlie preferential urn model 

In spite of the difhculty of occurrence of the conden- 
sation phenomenon in the Ehrenfest class with quenched 
disorder, it is possible to construct a model for the con- 
densation phenomenon. In order to research the conden- 
sation phenomenon in more detail, we consider a more 
concrete example, i.e., the disordered preferential urn 



model on a complete graph proposed in Ref. [19!]. The 
weight in the partition function is defined by 

p(/i„nO = (n,!)"'. (14) 

For this choice ofp{h,n), the condensation phenomenon 
can arise because 

oo 

H{h,z) = Y.^nl)^'''" (15) 

diverges when h = 1 and z > 1. 

The above choice of p{h, n) means the transition prob- 
ability u'{hi,ni) — {rii + 1)'*'. The above dynamics may 
not be adequate for physical systems, but this transition 
probability is related to the concept of preference, i.e., 
"the rich get richer" , which is important in the research 
field of complex networks [2^. Hence, it would be worth 
to consider such choice oip{h,n). 



B. Non-disordered case 

We here study the case with no disorder: All disorder 
parameters are hi = h = 1. Since the function H(h,z) 
of Eq. diverges with h = 1 and z > 1, we have to 
check the critical density. When h = 1, the saddle point 
equation ^ gives 



This means that the critical density pc in infinite. Hence, 
there is no condensation phenomenon, as described in 
Sec.III. 

It is easy to confirm that the condensation does not 
occur in the non-disordered case with h < 1. When h > 
1, a type of condensation could arise; we discuss the case 
in the following subsection. 



C. Disordered cases: Three types of condensation 

It has been revealed that the sites with the disorder 
parameter h.^ = 1 are important for the condensation 
phenomenon in the preferential urn model, but the non- 
disordered case is not sufficient for the occurrence of the 
condensation. Henceforth, we analytically and numeri- 
cally discuss three disordered cases which corresponds to 
different types of condensation phenomena, respectively. 



1. Complete condensation 

When there is a site hi > 1, Eq. (fT5)l has the radius of 
convergence (3 = 0, and the critical density Pc becomes 
zero from Eq. This condensation phenomenon shall 
be called the 'complete condensation'. In the complete 
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condensation, only one site contains almost of all parti- 
cles for any p, and the ratio of the number of particles 
containing the condensed site tends to be one in the limit 
of a large system. 

2. Genuine condensation induced by the disorder 

We here consider the disordered case in which each dis- 
order parameter h is assigned by the probability density 

= (l + a)(l -/i)" (0 < /i < 1). (17) 

When a = 0, i.e., we set 4'{h) as an uniform distribu- 
tion, we have no condensation phenomenon; the saddle 
point Zg gradually increases as p — > oo, but it converges 
to Zg = 1, so that the critical density is Pc = oo. 

Next, we check the case with a = 2 [23|. In this dis- 
ordered case, we get the critical density pc = 1.26, so 
that we have the condensation phenomenon in this dis- 
ordered case. We shall call the condensation as 'genuine 
condensation', and distinguish it from the complete con- 
densation and pseudo-condensation explained later. 

The occupation distribution in the steady state is ob- 
tained by Eq. . When we restrict the value of a in 
Eq. (fT7|) to an integer value and a > 0, we can calculate 
the approximate form of the occupation distribution as 

P{k) r^k-°'-^{\nk)-°'-^ (18) 

for large k region, by using Eqs. p2| . (fT4|. (fTT)) and set- 
ting Zs = 1. We note that the occupation distribution 
obtained analytically can not represent a condensed site; 
the occupation distribution is adequate except for the 
condensed site. When we have the condensation phe- 
nomenon, only a site deviates from the occupation dis- 
tribution obtained analytically. 

We here confirm the condensation phenomenon using 
numerical experiments. We have performed numerical 
experiments for the disordered urn model on a complete 
graph with lattice size L = 4000. All results are ob- 
tained by averaging 20 realizations. It is notable that 
the condensation arises only in the case in which there 
is at least one site with /i^ = 1, as discussed in Sec. III. 
Although we should assign the disorder parameter hi to 
each site by using the probability density 4'{h) defined 
as Eq. P?)) . a site with the disorder parameter hi = 1 
is difficult to emerge because we here treat a finite-size 
system. The value of hi near 1 is not sufficient to occur 
the condensation phenomenon, as discussed later. We 
therefore perform the following rescaling for the disorder 
parameters: 

1. Assign hi to each site by (t>{h) of Eq. P?|) . 

2. Search the maximum value of hi and set the value 

as /iniax ■ 

3. For all i, we perform the rescaling hi — > hi/h-an^^. 




k 

FIG. 2: The occupation distribution for the disordered case 
with a — 2. The solid line corresponds to the analytical result 
in Eq. (fT2)l with Za = 1. For details, see text. 



When we have the maximum value of /imax — 1; it is 
expected that the rescaling procedure has little influence 
on the disorder parameters {hi}. 

Figure [2] shows the results for a = 2.0. We have calcu- 
lated the cases with the density p — 1.4 and p = 2.4. The 
solid line in Fig. [5] corresponds to the analytical result of 
Eq. (fT2|) for Zg = 1. The occupation distribution shows 
the power law behavior in a certain k region, which is ex- 
pected from Eq. p^ . Note that because we performed 
the numerical experiments with L — 4000 and the results 
are averaged out over 20 simulations, the minimum value 
of P{k) is 1/(4000 X 20) = 1.25 x IQ-^. 

In Fig. [21 we can see that there are several points which 
deviate from the analytical result (solid line) because we 
take the average over 20 simulations. We here note that 
in each simulation, there is only one point which deviates 
from the analytical result. The deviating point means 
the fact that one site has quite a lot of particles in each 
simulation, which indicates the occurrence of the conden- 
sation phenomenon. Note that due to fluctuation among 
each simulation, the largest value of k in each simulation 
is different from each other, and hence we do not have a 
peak but the flat tail. 

We also see the condensation phenomenon by compar- 
ing with the case oi p = 1.4 and p = 2.4; the maximum 
value of k increases in the case of p = 2.4 compared with 
that of p = 1.4 in Fig. O However, the shapes of the 
distribution in small k region are the same in both cases. 
The analytical treatment suggests that the occupation 
distribution does not vary except for the condensed site 
in the condensed phase. Hence, we confirm the occur- 
rence of the condensation from the unchangeability of the 
occupation distribution in small k region and the shift to 
the right of the deviating points in the large k region. 

In order to see the condensation phenomenon more 
clearly, we have calculated the number of particles in the 
site which contains the maximum number of particles. 
The maximum number of particles in the condensed site 
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FIG. 3: The maximum number of particles in the site with 
condensation. The solid line represents Eq. (|19p . 



is expected to be 



N„ 



Lip- Pc), 



(19) 



because the system can contain only up to Lpc parti- 
cles without condensation. Figure [3] shows the numerical 
results for a — 2 and the analytical result of Eq. 
(the solid line). The maximum number iVmax shows a 
linear dependence on the density p. This result confirms 
that the macroscopic number (~ 0{L)) of particles is 
condensed only in a single site [2J| . 



3. Pseudo-condensation 




FIG. 4: The occupation distributions for the cases with 
single-site disorder. (a)The genuine condensation phe- 
nomenon (/idisordcr = 1.0). (b)Pseudo-condensation phe- 
nomenon (/idisorder = 0.8). The solid line corresponds to 
Eq. (f20)l . All data are averaged over 20 realizations. 



We here confirm that a site with the disorder parame- 
ter hi = 1 is necessary for the condensation phenomenon, 
as discussed in the previous subsection. We consider a 
single-site disordered case in which all disorder parame- 
ters hi = except for only one site with /idisordcr = 1- In 
this case, one can easily confirm that the critical density 
is equal to one (pc = 1) from Eq. The occupation 
distribution of Eq. ((T2)) becomes the Poisson distribution 

P{k)^e-P^^ (20) 

when the condensation phenomenon arises. We have 
performed numerical experiments with the lattice site 
L = 1000. 

Figure [IKa) shows the case in which only a single site 
has /idisordcr = 1- The numerical results are obtained 
by taking the average over 20 simulations. The solid 
line corresponds to the Poisson distribution of Eq. (PO)) . 
One can easily see the occurrence of the condensation 
phenomenon, as explained in Sec. IV. C. 2: compared with 
the case with p — 2.0, the maximum value of the occu- 
pation number increases in the case with p = 3.0. In 
addition, the occupation distributions are in good agree- 
ment with the analytical results, except for the condensed 
site. This condensation phenomenon corresponds to the 
genuine condensation, as described in Sec. IV. C. 2. 



In contrast, the genuine condensation phenomenon 
does not occur in the case where only a single site has 
/idisordcr = 0.8, as shown in Fig. |4l(b). One may see 
that the genuine condensation phenomenon occurs in 
Fig. SKb), because there are largely deviating points. 
However, in this case, the whole shape of the occupation 
distribution in small k region changes with the increase 
of the density p. When the genuine condensation arises, 
the increase of the density p does not change the occu- 
pation distribution in small k region, as explained in the 
previous subsection. The change of the occupation distri- 
bution means that we do not have the genuine condensa- 
tion phenomenon in spite of the existence of the deviating 
points in Fig. IH^b) . Evans et al. [1^ have called the sim- 
ilar type of condensation as 'pseudo-condensation'. In 
this case, the partition function is analytic even if the 
number of particles increases, so that there in no true 
phase transition. In summary, when there arc deviat- 
ing points in the occupation distribution and the whole 
shape of the occupation distribution varies with the in- 
crease of the density p, we call the phenomenon as the 
'pseudo-condensation'. It was revealed that the pseudo- 
condensation can occur in the preferential urn model de- 
pending on the nature of the quenched disorder. 
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V. CONCLUDING REMARKS 

In summary, the condensation phenomenon in 
'MaxweU-Boltzmann' statistics has been studied. We 
especially focused on disordered cases in the Ehrenfest 
class. We have developed the analytical treatment for 
the disordered cases, and clarified the conditions for the 
condensation. In addition, using the preferential urn 
model, we have demonstrated the condensation phenom- 
ena. It was revealed that only the disordered cases give 
the genuine condensation phenomenon in the preferential 
urn model. Furthermore, depending on the types of the 
quenched disorder, we have three types of condensation. 

In order to analyze the condensation phenomenon in 
disordered cases in more detail, we might need other an- 
alytical methods. For example, the canonical analysis 
proposed by Evans et al.^^ is useful for the study of 



the condensation in a finite size system; one can analyze 
the condensation transition and the structure of the con- 
densate, determining the precise shape and the size of 
the condensate in the condensed phase by the use of the 
canonical analysis. Although the grand canonical analy- 
sis for disordered cases has been developed in Ref. [13l |. 
we have not yet obtained an adequate analytical method 
for the disordered cases in the canonical ensemble. The 
disordered version of the canonical analysis would be im- 
portant for future research. 
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